Fault-tolerant quantum computation with high threshold in two dimensions 
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We present a scheme of fault-tolerant quantum computation for a local architecture in two spatial 
dimensions. The error threshold is 0.75% for each source in an error model with preparation, gate, 
storage and measurement errors. 
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Quantum computation is fragile. Exotic quantum 
states are created in the process, exhibiting entanglement 
among large numbers of particles across macroscopic dis- 
tances. In realistic physical systems, decoherence acts 
to transform these states into more classical ones, com- 
promising their computational power. Fortunately, the 
effects of decoherence can be counteracted by quantum 
error correction In fact, arbitrarily large quantum 
computations can be performed with arbitrary accuracy, 
provided the error level of the elementary components 
of the quantum computer is below a certain threshold. 
This is guaranteed by the threshold theorem for quan- 
tum computation 0, [3, Q, Hj] . 

Now that the threshold theorem has been established, 
it is important to devise methods for error correction 
which yield a high threshold, are robust against varia- 
tions of the error model, and can be implemented with 
small operational overhead. An additional desideratum 
is a simple architecture for the quantum computer, re- 
quiring no long-range interaction, for example. 

Recently, a threshold estimate of 3 x 10 -2 per opera- 
tion has been obtained for a method using post-selection 
jfj. An alternative scheme with high threshold com- 
bines topological quantum computation with state pu- 
rification 0. (See also [§].) In that approach, a sub- 
set of the universal gates are assumed to be error-free. 
Pure topological quantum computation ideally requires 
no error correction but often picks up a comparable poly- 
logarithmic overhead [9] in the Solovay-Kitaev construc- 
tion for approximating single- and two-qubit gates (c.f. 
(To|). fault tolerance is more difficult to achieve in archi- 
tectures where each qubit can only interact with other 
qubits in its immediate neighborhood. A fault tolerance 
threshold for a two-dimensional lattice of qubits with 
only local and nearest-neighbor gates is 1.9 x 10 -5 [111 ]. 

In this Letter, we present a scheme for fault-tolerant 
universal quantum computation on a two-dimensional 
lattice of qubits, requiring only a nearest-neighbor 
translation-invariant Ising interaction and single-qubit 
preparation and measurement. A fault tolerance thresh- 
old of 7.5 x 10 -3 for each error source is presented, with 
moderate resource scaling. This scheme is best suited 
for implementation with massive qubits where geometric 
constraints naturally play a role, such as cold atoms in 
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FIG. 1: (Color online.) The CNOT gate A(X) c , t (c: control, 
t: target) formed by topologically entangled lattice defects. 
Each pair of defects carries an encoded qubit. Defects exist as 
primal (blue) and dual (black) , and are created by local mea- 
surement. The primal correlation surface (light blue) shown 
here converts an incoming Pauli operator Z t into an outgoing 
Zt Z c , as required for a CNOT gate. 



optical lattices [12[ or two-dimensional ion traps [13[ . 

The presented scheme integrates methods of topologi- 
cal quantum computation, specifically the toric code [14[ , 
and magic state distillation [15[ into the one-way quan- 
tum computer (QCc) uM on cluster states. By employ- 
ing magic state distillation we improve the error thresh- 
old significantly beyond [13] , with the threshold value and 
overhead scaling now set by the topological error correc- 
tion. In this regard, we would like to emphasize that the 
three-dimensional cluster state is an intrinsically fault- 
tolerant substrate for quantum computation From 
the viewpoint of implementation it is desirable to reduce 
the spatial dimensionality of the scheme from three to 
two. To achieve this we turn the QCc into a sequential 
scheme in which the cluster state is created slice by slice. 

This Letter is organized as follows. First, we construct 
fault-tolerant universal gates for the QCc in three spa- 
tial dimensions. (See Fig. [1] for a CNOT gate.) Next, 
we perform the mapping to two dimensions. Finally, we 
present our error model and work out its threshold value. 

We consider a cluster state \4>)c on a lattice C with 
elementary cell as displayed in Fig. [2^. Qubits are located 
at the center of faces and edges of C. The lattice C is 
subdivided into three regions V, D and S. Each region 
has its purpose, shape and specific measurement basis 
for its qubits. The qubits in V are measured in the X- 
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FIG. 2: (Color online.) Lattice definitions, a) Elementary 
cell of the cluster lattice C. 1-chains of C (dashed lines), 
and graph edges (solid lines), b) A surface code obtained 
from a 2D cluster state by local X-measurements. c) A pair 
of electric ("e") or magnetic ("m") holes in the code plane 

each support an encoded qubit. z e ^ m and X e ^ m denote the 
encoded Pauli operators Z and X, respectively. 



basis, the qubits in D in the Z-basis, and the qubits in 
S in either the F-basis or the eigenbasis (X -\-Y)/ y/2. 

V fills up most of the cluster. D is composed of thick 
line-like structures, named defects. S is composed of well- 
separated qubit locations interspersed among the defects. 
As described in greater detail below, the cluster region 

V provides topological error correction, while regions D 
and S specify the Clifford and non-Clifford parts of a 
quantum algorithm, respectively. 

We can break up this measurement pattern into gate 
simulations by establishing the following correspondence: 
quantum gates <-> quantum correlations <-> surfaces, as il- 
lustrated for the CNOT gate in Fig. CD The first part of 
this correspondence has been established in [16]. For 
the second part homology comes into play. The cor- 
relations of \<j))c (i- e -? the stabilizers) can be identified 
with 2-chains (surfaces) in £, while errors map to 1- 
chains (lines). Homological equivalence of the chains im- 
plies physical equivalence of the corresponding operators 
[171 ]. This correspondence is key to the presented scheme. 
Gates are specified by a set of surfaces with input and 
output boundaries, and syndrome measurements corre- 
spond to closed surfaces (having no boundary). 

Formally, C is regarded as a chain complex, C = 
{C 3 ,C 2 ,Ci,C }. It has a dual C = {C 3 ,C 2 , Ci,C } 
whose cubes q G C3 map to sites Co G Co of £, whose 
faces c 2 G C 2 map to edges c\ G C\ of £, etc. The chains 
have coefficients in Z 2 . One may switch back and forth 
between C and C by a duality transformation *( ). £, C 
are equipped with a boundary map <9, where d o d = 0. 

Operators may be associated with chains as follows. 
Suppose that for each qubit location a in a chain c, 
a G {c}, there exists an operator £ a , with [£ a ,£&] = 
for all a, 6 G {c}. Then, we define E(c) := E[aG{c} 
Cluster state correlations are associated with 2-chains. 
Specifically, all elements in the cluster state stabilizer 
take the form K(c 2 )K(c 2 ) with c 2 G C 2 , c 2 G C 2 , and 
K(c 2 ) = X(c 2 )Z(dc 2 ), K(c 2 ) = X(c 2 )Z(dc 2 ). Only 
those stabilizer elements compatible with the local mea- 



surement scheme are useful for information processing. 
In particular, they need to commute with the measure- 
ments in V and D, 

[K(c 2 )K(c 2 ),X a } = 0, a eV, m 
[K(c 2 )K(c 2 ),Z b ] =0, beD. [l> 

This condition may again be expressed in terms of the 
chains c 2 , c 2 directly, which we will do below. 

Topological error correction in V . Inside V the con- 
straint (pQ) implies dc 2 = 0, dc 2 =0. In particular, 
these conditions are obeyed for c 2 = dcs, c 2 = dcs. 
For each elementary cube q G C3, q G C3 the cluster 
stabilizers K(dq), K(dq) can be measured by the local 
X-measurement and classical post-processing. 

The optimal error correction procedure for V can be 
mapped to a model from classical statistical mechanics, 
the random plaquette 7L 2 - gauge model in three dimensions 
(3D-RPGM) [18], for which a fault tolerance threshold 
of 3.3 x 10 -2 for local noise has been found in numer- 
ical simulations 0. (See also [20|.) Here we use the 
minimum weight chain matching algorithm [2l| for er- 
ror correction. It yields a slightly smaller threshold of 
2.9 x 10 -2 [22[ but is computationally efficient. Various 
error sources eat away at this 3% error budget. 

Cluster states and surface codes. The connection be- 
tween a 2D cluster state and a surface code is illustrated 
in Fig. [2)3. The extra spatial dimension in a 3D cluster 
state allows to evolve coded states in "simulated time". 
The number of qubits which can be encoded in a surface 
code depends solely on the surface topology. Here we 
consider a plane with pairs of either electric or magnetic 
holes; see Fig.[2t. A magnetic hole is a plaquette / where 
the associated stabilizer generator S\j(f) = Z(df) is not 
enforced on the code space, and an electric hole is a site 
s where the associated stabilizer S+(s) = X(d #s) is not 
enforced on the code space, where denotes the dual- 
ity transformation in 2D. Each hole is the intersection of 
a defect strand with a constant-time slice. 

A pair of holes supports a qubit. For a pair of magnetic 
holes /, the encoded spin flip operator is X = X(ci), 
with {dci} = {^f,^f'} : and the encoded phase flip op- 
erator is Z m = Z(ci), with ci = df or c x ^ df. The 
operator Z(df + df f ) is in the code stabilizer. For a pair 
of electric holes 5, s f we have X = X(c[), with c[ = d#s, 
Z e = Z(ci), with {da} = {s, s 7 }, and X(d*s + d*s') is 
in the code stabilizer. 

Quantum logic. The CNOT gate is realized by link- 
ing primal and dual defects as displayed in Fig. [TJ To 
explain the functioning of the gate we refer to Theorem 
1 of [l6[ • We consider a block shaped cluster C where the 
elementary cell of Fig. [2^ is repeated an integer number 
of times along each direction. One of these directions is 
singled out as "simulated time". The two perpendicu- 
lar slices of the cluster at the earliest and latest times 
contain the supports / and O for the encoded input and 
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FIG. 3: (Color online.) Remaining gates for universal fault- 
tolerant computation. The relevant correlation surfaces are 
shown in light blue and gray. Replace Out (put) by In(put) 
for a measurement, a) Preparation of a Z-eigenstate for an 
electric qubit. b) Preparation of an X-eigenstate for an elec- 
tric qubit. c) Creation of a Bell pair among a bare S'-qubit 
and an encoded qubit. 
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FIG. 4: Elementary cell of the 2D lattice. Temporal order 
of operations in V: The labels on the edges denote the time 
steps at which the corresponding A(Z) gate is performed. The 
labels at the syndrome vertices ("o") denote measurement 
and (re-) preparation times [£m, £p], and the labels at the code 
vertices ("•") denote times for Hadamard gates (tH,t' H ). The 
pattern is periodic in space, and in time with period six. 



output qubits, respectively, with I,Oc {C\} encoded by 
the surface code of Fig. [2k. 

The set M on which the measurement pattern is de- 
fined (c.f. Thm 1 of [l6j]) is composed of V and D, 
M = V U D. Due to the presence of a primal lattice C 
and a dual lattice £, it is convenient to subdivide the 
sets V and D into primal and dual subsets. Specifi- 
cally, V = V P U Vd, with V p C {C 2 }, Vd C {C 2 }, and 
D = D P U D d , with D p C {Ci}, D d C {d}. 

With these definitions, we can now prove the func- 
tioning of the CNOT gate in Fig. [H The gate cluster C 
contains the regions V p , Vd-, D p , Dd, I and O. In this 
setting, condition (pQ) implies for the correlation surfaces: 

{c 2 } C V p , {dc 2 } C DpUlUO, ( . 

{c 2 } C VdU/UO, {dc 2 } C D d . 1 } 

One such (primal) correlation surface is depicted in 
Fig. [TJ The corresponding stabilizer of \<j))c, after mea- 
surement of the qubits in M = VUD, implies a stabilizer 
±XjX,oX,o = ±ZliKX)cX,oKX)l t for |*) jo . 
Three similar surfaces imply the stabilizer elements 

±X,iX,o, ±XjX,o and ±XjX,oX,o for |*>/o- 
Theorem 1 of [16[ is applied with U = A(X) Cyt . □ 
Further elements of a fault-tolerant QCc-computation 
are shown in Fig. [3j Fault-tolerant preparation of en- 
coded X- and Z-eigenstates for the electric qubits are 
displayed in Figs. [3^t and [3)3, which can be reversed to 
denote measurements. These operations, together with 
the CNOT gate of Fig. [U comprise the set of topologi- 
cally protected gates. Fig. [3b shows the creation of a Bell 
pair between a bare S'-qubit and a qubit encoded with a 
surface code (electric). The shown correlation surfaces 
c 2 , c 2 are such that {02} C V p , {dc2} C D p U S U O, 
{02} C Vd U S U O, {dc2} — 0. The corresponding sta- 
bilizers K(c2), K(c2) imply, after local measurement of 
the qubits in V and D, the stabilizer generators ±ZsZo, 
±X S X for the state \^)so- Thus, \SP)so is a Bell state 
with the qubit located on O being encoded. Measure- 
ment of the bare qubit on S in the eigenbasis of Y or 



(X + Y)/>/2 yields on O an encoded state \Y) = |0)+i|l) 
or \A) = |0) + e m / 4 |l), respectively. These states are 
noisy and therefore subsequently purified via magic state 
distillation [15[. Finally, they are used in teleportation 
circuits (see Fig. 10.25 of |24() to generate the fault- 
tolerant gates exp(z7r/4X) and exp(i7r/8Z). This com- 
pletes the universal fault-tolerant gate set. 

Mapping to the 2D lattice. The dimensionality of the 
spatial layout can be reduced by one if the cluster is 
created slice by slice. That is, we convert the axis of 
"simulated time" — introduced as a means to explain the 
connection with surface codes — into real time. 

Cluster qubits located on time-like edges of C or C be- 
come syndrome qubits, which are periodically measured. 
Qubits on space-like edges become code qubits. Time- 
like oriented A(Z) gates are mapped to Hadamard gates, 
while space-like oriented A(Z) gates remain unchanged. 

The temporal order of operations is displayed in Fig.[H 
Note that every qubit is acted upon by an operation in 
every time step. The mapping to the two-dimensional 
structure has no impact on information processing. In 
particular, the error correction procedure is still the same 
as in fault-tolerant quantum memory with the toric code. 

Error model and threshold. There are two separate 
thresholds, one for the Clifford operations and one for 
the non-Clifford operations. The former threshold de- 
rives from topological error correction and the latter from 
magic state distillation. The overall threshold is set by 
the smaller of the two. 

Mapping to a single-layer 2D structure slightly mod- 
ifies the effective error model on the lattices C and £, 
as compared to [17]. Specifically, we assume the fol- 
lowing: 1) Erroneous operations are modeled by perfect 
operations preceded or followed by a partially depolariz- 
ing single- or two-qubit error channel T\ = (1 — + 
pi/3 ([X] + [Y] + [Z]), T 2 = (l-p 2 )[I]+p 2 /l5 ([X a X b ]+..+ 
[Z a Zb]). The error sources are a) the preparation of the 
individual qubit states |+) (error probability pp), b) the 
Hadamard gates (error probability pi), c) the A(Z) gates 
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(error probability P2), d) measurement (error probability 
Pm)- 2) Classical syndrome processing is instantaneous. 

When calculating a threshold, we assume that all er- 
ror sources are equally strong, p\ = P2 — Pm — Vp '•— P- 
Storage errors need not be considered because no qubit 
is ever idle between preparation and measurement. This 
model encompasses realistic error sources such as local in- 
homogeneity of electric and magnetic fields, fluctuations 
in laser intensity, and imperfect photodetectors. 

The topological threshold for each physical source is 
estimated by numerical simulations to be 

p c = 7.5 x 1(T 3 . (3) 

A similar threshold persists under modifications of the 
error model such as higher weight errors 

Regarding the distillation threshold, the residual error 
€1 at level I undergoes the recursion e/ — > q + i = 35q 3 
(to leading order) [151 ]. The initial distillation error eo 
arises through the effective error on an .S-qubit, with 
eo = 6p. The distillation threshold p c for each physi- 
cal error source is then p c = 1/6^35 ~ 2.8 x 10 -2 . The 
purification threshold is much larger than the topological 
threshold, and therefore the overall threshold for fault- 
tolerant QCc-computation is given by Eq. ([3|). 

Overhead, fault tolerance leads to a poly-logarithmic 
increase of operational resources. Both the overheads in 
topological error correction and in magic state distillation 
are described by a characteristic exponent: j t op — 3 and 
7 ms = log 3 15. The larger one dominates the resource 
scaling. Given bare circuit size S, the encoded circuit 
size S' scales as S' ~S log 3 S. 

Conclusion. We have presented a scheme of fault- 
tolerant quantum computation in a two-dimensional lo- 
cal architecture with high error threshold and mod- 
erate overhead in resource scaling. The threshold of 
7.5 x 10 -3 is the highest known for a local architec- 
ture. Our scheme only requires local and translation- 
invariant nearest-neighbor interaction in a single-layer 
two-dimensional lattice. Small-scale experimental de- 
vices may be realized in optical lattices, segmented ion 
traps, or arrays of quantum dots or superconducting 
qubits where short-range interaction is preferred. 
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